CERN-TH/2001-286 
DESY 01-163 
IFIC/01-57, FTUV/011019 

1 
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We compute meson propagators in finite- volume quenched chiral perturbation theory. 



Recent progress in the formulation of chiral 
fermions on the lattice makes it possible to ap- 
proach the regime of the physical light quark 
masses. However, lattice volumes V = L'^ needed 
to ensure negligible finite size effects are pro- 
hibitively large for such light quarks. Only when 
L is very large compared to the Compton wave- 
length of the lightest particles of the theory (i.e. 
the pions) are the volume effects exponentially 
suppressed. Fortunately chiral perturbation the- 
ory (xPT) can predict the volume effects that oc- 
cur in the regime m^L < 1. In this way the fi- 
nite volume becomes a distinct advantage, and a 
finite-size scaling analysis is actually a very useful 
tool to extract the physical low energy constants. 
This has been demonstrated in simpler models 
0. 

The xPT predictions in full QCD for the fi- 
nite size scaling of quantities such as the chiral 
condensate or the propagators of scalar and vec- 
tor densities have been computed in |^,^. Un- 
fortunately, lattice simulations in the regime of 
very light quark masses are presently only possi- 
ble in the quenched approximation, and the pre- 
dictions of chiral perturbation theory are quite 
drastically modified in this approximation [Q . We 
present here some results for the propagators of 
the scalar and pseudoscalar densities at finite vol- 
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ume in quenched xPT (QxPT). Full details will 
be presented in a forthcoming publication. 

The low-energy limit of QCD with light quarks 
is a chiral Lagrangian of the Goldstone bosons re- 
sulting from spontaneous chiral symmetry break- 
ing. In an infinite volume, the effective La- 
grangian is approximated by an expansion in 
powers of the pion momentum p and mass m^. 
This is the standard p- expansion. In a finite vol- 
ume the p-expansion is still good if the Compton 
wavelength of the particle is much smaller than 
the size of the box ( L ^ ^ ). In the oppo- 
site limit, L ^ the p-expansion breaks down 
due to propagation of pions with zero momenta 
1^ . A convenient expansion for this regime is the 
so-called e- expansion in which rrin ^ p^ ^ 
1^. By the Gell-Mann-Oakes-Renner relation 
this corresponds to keeping pii = rriiTiV of or- 
der unity while V is taken large. The zero 
mode of the pion can be isolated by factorizing 
U{x) = Uq exp i^/2^{x) / into the constant col- 
lective field Uo and the pion fluctuations £,{x). 
The difficulty then comes from the fact that the 
integral over Uq needs to be done exactly, while 
ordinary xPT applies to the non-zero mode inte- 
gration. At leading order one obtains the suitably 
normalized partition function 



Z = I dUo exp 

lsU{Nf) 



-VTtM {Uo + U+ 
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where M is the quark mass matrix. It is interest- 
ing to consider averages in sectors of fixed topol- 
ogy v as well B . To the same order 



Zy~ i rfC/o (dctJJo)" exp 

JU{Nj) 



-VTtM (Uo + C/n* 



The e-expansion can be worked out also in the 
case of QxPT. We have considered two methods 
to take the quenched limit. The first is the su- 
persymmetric method suitably refined to be 
valid at the non-perturbative level Here the 
fermionic determinant of iV„ valence quarks is 
cancelled by adding to the theory N^^ flavors of 
ghost bosonic quarks. Assuming a supergroup 
generalization of the chiral symmetry breaking 
pattern of QCD, the field U{x) becomes an el- 
ement of a graded group with both Goldstone 
bosons and Goldstone fermions. The integration 
over these fields is taken on what is called the 
maximal Riemannian submanifold of the super- 
group Gl(A^i,|A'^„) This is a combination of 
compact (as for the usual Goldstone bosons) and 
non-compact integration domain (for the Gold- 
stone bosons of the ghost-ghost block). The fac- 
torization of U{x) is still as above, but ^(x) is no 
longer traceless since the singlet field does not de- 
couple in the quenched approximation. To lead- 
ing order in the e-expansion the partition function 
in a fixed topological sector is 



= JdUod^ (SdetC/o)''exp|-FStrM([/o + C/, 
id^x^ 



can only be performed through series expansions 
As a first check, we have confirmed the cal- 
culation of the leading correction to the chiral 
condensate in the e-expansion now using the 
supersymmetric method. Likewise, as a check on 
the results presented below we have performed 
all calculations both ways, and thus compared 
the resulting series expansions from the replica 
method with the closed expressions obtained by 
the supersymmetric method. In all cases we find 
complete agreement. 

Here we present the first quenched results for 
correlation functions in the e-expansion. We be- 
gin with the scalar and pseudoscalar correlation 
functions s"(a;) = (5'°(a;)5"(0)) and p°(a;) = 
(P'^(a;)F'^(0)) to ©(e^), which is the leading order 
contribution to the space-time dependent terms. 
The 0(1) contributions are constant. In a sector 
with fixed topology we find for the space integrals 
of flavor singlets 



s"(t) = const - — r G(t)a- 
p%t) = const + — G(^)a^ 



where V = L^T and 

with the chiral condensate |q 



a_ - 4 



with ^{x) 
trace. 



Str^(x) and Str denotes the super- 



^=/i(/.(/i)i^,(A*) + /.+i(A*)if.-i(M)) + -. 



An alternative to the supersymmetric method 
is the replica method , which amounts to taking 
the Nf ^0 hmit of the full QCD result. Chiral 
perturbation theory done this way is equivalent 
to that of the supersymmetric method, but the 
Feynman rules are somewhat simpler. Its short- 
coming in the present context is that the inte- 
gral over the zero momentum modes in general 



The functions G{t) and A{t) are, with r = ^, 



A{t) 



G{t) = -^T'h,{r) + ^Th,{r) 



hi{T) 
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h2{T) = — 

^ ' 24 



^ ^ 30 



The origin of the functions /ii.2(t) has been dis- 
cussed in detail in the hterature [|0[. As a check 
on the above resuh, in the hmit /i — > oo wc find 



2S2 



m)-A{t)] 



which coincides with the leading order result from 
the p-expansion. This is simply tree level propa- 
gation of the flavor singlet. 

We have also computed the flavored correla- 
tion functions. The calculation is more involved 
because in the supersymmetric formulation it is 
necessary to work with at least Ny = 2 since 
the sources contain two flavors. All details will be 
given in a forthcoming publication. Here we just 
show a plot of the pseudoscalar correlation func- 
tion p°'(t), a = 1, 2, 3 in the e and p-expansions, 
for different values of the quark mass close to 
the range of validity of both regimes^. As in 
the unquenched case, there is an overlap region 
at MjrL ~ 1 where both expansions are good if 
mEt/ > 1. 

The volume dependence of the above formu- 
lae provides an extra handle to extract the low- 
energy constants S, i^jri ™o ^^'^ ^ i'^ ^^^^ regime. 
This technique has already been successfully ap- 
plied to extract the constant E from the spectral 
density of the Dirac operator and from the quark 
condensate at finite volume |12l . 
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Figure 1. 
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Flavor- non-singlet p°(t) in lattice 
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Shown is the p-expansion (dashed lines, m = 
10, 3, 0.48MeV) and the e-expansion (solid lines, 
for m = 0.48, 0.2MeV). As m is decreased the 
curves go upward, as expected. In the intermedi- 
ate mass region (m ~ 0.5Mev) there is nice agree- 
ment between the two expansions. 
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